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■ Abstract. We compare parameters of well-observed type II SN1999em de- 
' rived by M.Hamuy and D.Nadyozhin based on Litvinova & Nadyozhin (1985) 

^ |. analytic fits with those found from our simulations using our radiative hydro 

D ' code STELLA. The same code applied to models of SN1993J allows us to es- 

[ timate systematic errors of extracting foreground extinction toward SN1993J 

suggested by Clocchiatti et al. (1995) which is based on the assumption of black 
body radiation of the supernova envelope near the first maximum light after 
I shock break out. A new implicit two-temperature hydro code code SUPREMNA 

^-H . is introduced which self-consistently takes into account the kinetics of ioniza- 

^ ' tion, electron thermal conduction, and radiative losses. Finally, a combination 

^sD i of STELLA and SUPREMNA allows US to use the same SNIa models both for build- 

0^ ■ ing their light curves and predicting X-ray spectra of young Supernova remnants 

CO ' such as Tycho and Kepler. For the comparison of theoretical results with the 

0^ I observations we used data on Tycho SNR obtained with XMM~Newton space 

. telescope. 
. 

o ■ 

^1 1. Introduction 

H ■ 

I Three topics will be discussed briefly here. First, we show a way for extracting 

■ SN II parameters from UBV light curves (using SN1999em as an example). Then 
^ I we turn to estimating foreground extinction (example of SN1993J). And finally, 

• ■ we demonstrate using SN la models for predictions both of UBVRI light curves 

rS ' of SNe and X-ray spectra of young SNRs. 

2. SN II light curve theory vs. observations 

Nadyozhin (2003) and Hamuy (2003) have obtained masses M, radii R, and 
explosion energies E for a set of SNe II. Both have taken the observed values 
of My, At, Mph on the light curve plateau and extracted A4, R, E, based on the 
relations found by Litvinova & Nadyozhin (1983, 1985), — LN85 hereafter. 

Let us look into parameters found by Nadyozhin (2003) and Hamuy (2003) 
for an example of well-studied SN 1999em. Note that Nadyozhin's (2003) Eq.(l) 
changed units of E from LN85 and 2 typos are introduced by MN into this paper; 
see the original version astro-ph/0303411. 

From the Table we see that the discrepancy of derived parameters for the 
same supernova is large. While input parameters differ within 10% only, the 
derived ones may be in disagreement by a factor of 2. There is an appreciable 
difference in the assumed distance D in the two papers, and hence in My. The 

1 



2 



Blinnikov et al. 



Input 


iNaclyoznin [Zyjvoj 


xiamuy [Zuuo) 




D 


12.38 


10.7 


Mpc 


Mv 


-16.78 


-16.44 




At 


110 


124 


days 




2900 


3290 


km/s 


Derived 
E 
M 
R 


0.63 
13.2 

569 


1 9-1-0.6 
-^■^-0.3 


foe 

Mq 
Rq 



values of E, A4, and R scale with the distance as: 

Thus, it is very important to know D with as high accuracy as possible. But the 
difference may be explained by different D only for R, but not for E and Ai: 
other small differences are also important. Hence, the LN85 process is rather 
unstable. 

If one looks into the UBV light curves in Litvinova & Nadyozhin (1983) 
one can see that they are all similar with a pronounced plateau, because they 
were built in a simple equilibrium diffusion approximation, while well observed 
SNe II have a good plateau only in V and redder filters, but not in B and U - 
see, e.g. Leonard et al. (2002) and Fig^ below. 

This means that more detailed models are needed in order to reproduce the 
observed UBV behavior. We have considered an example of the well-studied 
SN 1999em with a goal to extract its parameters using more sophisticated mod- 
ern multi-group SN II LC modeling. We have used our code Stella (Blinnikov 
et al. 1998; Sorokina & Bhnnikov 2003, 2004). See Baklanov (2002) where a 
theoretical catalog of UBV, M^oi light curves is built for very different SN II 
presupernovae for a set of A4,R,E. In addition to those parameters used by 
LN85 the mass A^Ni of radioactive ^^Ni is also accounted for. 

To model SN 1999em we construct non-evolutionary preSN models with 
extended hydrogen-rich envelopes, similar to LN85, but a bit more realistic 
because they have a compact massive core and ^^Ni mixed to the envelope as 
in the models for SN 1987A (Blinnikov et al. 2000). The envelope has mass 
fraction of H, X=0.7, and the fraction of metals was assumed Z=0.03 for the 
first models (in proportion to cosmic distribution with this Z). 

The Fig^ shows our model light curves when the presupernova parameters 
(and reddening) are those of Hamuy (2003) & Nadyozhin (2003). We see that 
fluxes on the plateau disagree up to two stellar magnitudes with observations. 

In search of best-fitting model were have varied the input parameters in wide 
range and found a reasonably good fit for M = 15Mq, R = 450Rq, E = 7x lO^'' 
ergs. Mm = OMMq (FigH left). This set is good only for D = 7.5 Mpc (EPM, 
Hamuy et al. 2001), which is much shorter scale than the Cepheid distance 
found by Leonard et al. (2003). A very nice fit is found for the same parameters, 
but for low metal abundance, Z=0.004 (FigEl right). Not only UBV, but also 
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Figure 1. UBV for SN 1999em parameters — Hamuy (left), Nadyozhin (right) 
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Figure 2. UBV short distance scale: D = 7.5 Mpc, metals Z^^O.OS (left), 
Z=0.004 (right). 



the dependence of photospheric speed Uph{t) is good for the last model. Low 
metal abundance agrees well with the results of spectral modeling by Baron et 
al. (2000). For the long distance scale one can get a satisfactory UBV only for 
huge R ~ 10^i?Q. Work on this is in progress. 



3. Extinction toward SN 1993J 

The determination of foreground extinction toward supernovae (SNe) is very 
important. A technique for estimating the extinction has been suggested by 
Clocchiatti et al. (1995) which is based on the assumption of black body SN 
radiation near the first maximum light after shock breakout. Here we note 
that the true spectrum of an SN after the shock breakout, theoretically, is not 
a pure black body (due to importance of scattering and non-gray absorption 
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in the ejecta). Using spectra which were calculated by the Stella code for 
the 13C model of SN 1993J (Woosley et al. 1994; Bhnnikov et al. 1998) and 
which fit energy distribution in the observer frame sufficiently well, we found 
that the values of Ay may differ substantially from found under the black body 
assumption (Fig|21). 




log(L'[Hz]) lime, days 

Figure 3. Left. 1.12 days after SN explosion. Solid: model. Dashed: black 
body fitted to the model in the visible range (labeled 'FIT') with Cardelli et al. 
(1989) absorption law. The best fit gives T = 39300 K and Ay = 0.241 (while 
true absorption in this numerical experiment is zero!). Right. Dependence of 
faked extinction parameter Ay on time for different extinction laws: solid for 
a simple exponential law and dashed for CardeUi et al. (1989) 



4. SN la models: light curves and Tycho Supernova Remnant 

Here we concentrate on 3D SN la models computed at MPA (Reinecke, Hille- 
brandt, & Niemeyer 2002). Our assumptions on their light curve models using 
STELLA are described in Sorokina & Blinnikov (2003), Blinnikov &: Sorokina 
(2004). The composition of early MPA models was rather crude. Recently 
Travaglio et al. (2004) computed a detailed explosive nucleosynthesis for some 
of them. We show Flight curves for those models and a bolometric flux for cen- 
trally ignited c3_3d_256 in Fig0J The curves are in reasonably good agreement 
with observations. 

It is interesting to compare the properties of young SNRs, such as Tycho, 
with predictions of the same models that have been used for LC modeling. For 
the first hundreds years the pattern of flow has a forward shock sweeping up ISM 
material and a reverse shock propagating into the slower-moving ejecta which 
we call Nadyozhin-Chevalier stage: a self-similar solution by Nadyozhin (1981, 
1985) and Chevalier (1982) is much more realistic than widely used Sedov solu- 
tion. For this stage Sorokina et al. (2004) have developed a new ID Newtonian 
implicit hydro code supremna which allows us to study dynamics and radiation 
of young SNRs without using approximate solutions. 
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Figure 4. Left. V light curves for new MPA models. Right. Bolometric 
light curves for c3_3d_256 (dotted Lboi, solid Luvom) 



The main features of the code: ionization energy is taken into account in 
the equation of state; electron thermal conduction with saturation at the sound 
speed and radiative losses are not neglected; Te 7^ Tj (equilibration by Coulomb 
collisions). The X-ray spectrum is computed by a time-dependent ionization 
code which calculates evolution of the ionization stages for every mesh zone at 
each hydro time step for all ions of 15 elements; it takes into account collisional 
ionization, autoionization, photorecombination, dielectronic recombination, and 
charge transfer. Comparison with other codes used for young SNRs is given in 
the table. 
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As an example, we show in Fig 


3 the X-ray spectra for the age 



430 years for W7 model and for the same MPA model MRO and for the same 
parameters of ISM as used by Blinnikov & Sorokina (2004) and Kosenko et al. 
(2004). The ID W7 model has a problem reproducing iron emission without 
an additional mixing, as noted already by Itoh et al. (1988), while MPA models 
are well mixed from the beginning and they show a prominent iron line. There 
may be a problem with lines of intermediate mass metals, such as Ca and Si, 
but the development of thermonuclear 3D deflagration models may resolve this 
problem. 

An important result of our work is the development of thermal instability in 
metal-rich ejecta in the models with low electron conduction. This confirms the 
old result found by Hamilton & Sarazin (1984) which was ignored in numerical 
models of young SNRs. 
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Figure 5. X-ray spectrum for W7, convolved with the response matrix of 
EPIC PN on XMM-Newton (left), the same for MPA model MRO (right). 



5. Conclusions 



• To get SN parameters, one has to craft a full light curve model and Vph{t). 

• In case of SN 1999em we easily find reasonable parameters and low Z for 
EPM distance of 7.5 Mpc. The Cepheid distance requires a larger R of 
the presupernova. 

• Using black body assumption for early SN spectra may introduce a sys- 
tematic error into estimates of Ay- 

• New MPA 3D models of SN la produce faster bolometric LCs, closer to 
observed ones. 



• Current 3D deflagration MPA models can easily explain iron in spectra of 
young SNRs. 

• Energy losses may be very important in the evolution of metal-rich ejecta 
of SN la. They must be taken into account in realistic simulations of their 
dynamics. 
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